Evolution of a localized electron spin in a nuclear spin environment 
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Motivated by recent interest in the role of the hyperfine interaction in quantum dots we study the 
dynamics of a localized electron spin coupled to many nuclei. An important feature of the model 
is that the coupling to an individual nuclear spin depends on its position in the quantum dot. We 
introduce a semi-classical description of the system valid in the limit of a large number of nuclei 
and analyze the resulting classical dynamics. Contrary to a natural assumption, the correlation 
functions of electron spin with an arbitrary initial condition show no decay in time. Rather, they 
exhibit complicated undamped oscillations. This may be attributed to the fact that the system 
has many integrals of motion and is close to an integrable one. The ensemble averaged correlation 
functions do exhibit a slow decay (~ l/ln(t)) for t — > oo. 

PACS numbers: 73.21.La,85.35.Be,72.25.Rb 



I. INTRODUCTION 

The coherent manipulation of localized spins in solid 
state systems is currently a very active field research 1 . 
One of the most ambitious goals in this field of re- 
search is developing a quantum bit, or qubit. Such 
qubits would form the basic building blocks of quantum 
computers2i2i£i£. The strength of the qubit coupling to 
the environment should be sufficiently well controlled so 
that the residual environmental interactions are below a 
threshold set by the operating speed of the devic o 6 i 7 i 8 . 

Since the first proposal in Ref. \% the qubit candidate 
based on localized electron spins in quantum dots has 
been the subject of many theoretical and experimental 
studies. There have already been quite a few studies 
of the relaxation time of electron spins in GaAs quan- 
tum dots. The role of various spin-orbit related mech- 
anisms has been considered for spin-flip transitions be- 
tween singlet and triplets^ and Zeeman split doublets^. 
The hyperfine mediated spin-flip rates were investigated 
for singlet-tripleliii and doubled transitions. A novel 
spin-flip mechanism due to interface motion has also been 
investigated^ and Coulomb blockade effects in the nu- 
clear spin relaxation time were looked at in Ref. IhH 
For GaAs quantum dots measurements of spin relax- 
ation times, or T%, between singlet-triplet statesi^*^ and 
doubled states have been done. The measurements give 
Ti w 200 ^s (singlet-triplet) and T x > 50 ^s (doublet). 
For the doublet measurements only a lower bound for 
the relaxation time is obtained^, so comparison to vari- 
ous theoretical calculations remains difficult. 

The operation of a qubit requires a coherent superpo- 
sition of states that persists during the time of the qubit 
operation. Several papers have addressed the question 
of the decoherence time of electron spin in GaAs quan- 
tum dots. The decoherence of a single electron spin may 
be caused by inhomogeneous hyperfine interactioniSiiS 
but when an ensemble of electron spins is considered the 
timescale of decoherence is determined by the fluctua- 
tions in the total field which acts on the electron spin 
due to the nucleiSSiSi. The role of fluctuations in the 



nuclear spin system, socalled spectral diffusion, on the 
electron spin decoherence has also been studied^ 2 *^. 

In Ref. the characteristic timescale for the decay of a 
specific correlation function was associated with the deco- 
herence time of the electron spin. The decohence in this 
case was attributed to the spatially dependent hyperfine 
coupling constant which caused small frequency changes 
through flip-flop processes involving spatially separated 
nuclei. A different approach was used in Ref. |2(| which 
relied on representing the nuclear spin system as an ef- 
fective nuclear magnetic field that couples to the electron 
spin in the quantum dot. Merkulov et al. discussed some 
basic features of this semi-classical approach but the cal- 
culations were done for ensemble averaged quantities 20 . 

In the present work we combine the approaches 
mentioned and extend the semiclassical effective-field 
methodic to include the effects of the spatially varying 
hyperfine coupling constant. Due to the large difference 
of timescales for the electron and nuclear spin systems 
we are able to solve the problem in two steps. In the first 
step we establish that the nuclear system can be treated 
as an adiabatic effective nuclear magnetic field acting on 
the electron. The latter step involves the back action of 
the electron spin which will determine the evolution of 
the nuclear spins. The fact that the single electron spin 
is coupled to a large number of nuclei in the quantum 
dot, N, but each nucleus is only coupled to the single 
electron spin hints at an asymmetry in the behavior of 
the electron and nuclear spins. The electron precesses 
in an effective nuclear magnetic field which is due to the 
whole nuclear spin system. This field is oc y/W times 
larger than the corresponding hyperfine field due to the 
single electron spin in which the nuclear spins precess. 
Thus the dynamics of the electron are much faster than 
the dynamics of the nuclei. Also, the large number of 
nuclei involved makes it possible to treat the nuclear sys- 
tem in a semi-classical waj*i2iS2*24. We do not include the 
dipole-dipole interaction of the nuclei, since the effects of 
such interaction can only become noticeable at very long 
timescale ( idip — 10 _3 sec, for GaAs quantum dot o 20 i 25 ). 

The resulting dynamics are represented by a set of 
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equations of motion for several subsystems of nuclear 
spins. Each subsystem is characterized by the same value 
of coupling to the electron spin. The dynamical equations 
are non-linear and comprise many degrees of freedom. 
From this, one generally expects chaotic and ergodic dy- 
namics, so that the memory about initial conditions is 
lost at a certain time scale. This would result in the de- 
cay of correlation functions at this time scale and was 
an initial motivation of this research. However, we prove 
that the actual dynamics is not chaotic. The calculated 
correlation functions show complex, yet regular, oscilla- 
tions that persist for long timescale without discernible 
decay. We explain this by conjecturing that the system 
has many integrals of motion so it is close to exactly in- 
tegrable one. 

If we average the dynamics over all possible initial con- 
ditions, the averaged correlation function does show de- 
cay in time. For large times, it is inversely proportional 
to logarithm of time. 

The remaining text of the paper is organized as fol- 
lows. In Sec. [HI the hyperfine interaction in a quantum 
dot is presented and its semi-classical representation dis- 
cussed. The adiabatic approximation for the electron 
spin dynamics is discussed in Sec. IIIII and some analytic 
properties of the dynamical equations are presented in 
Sec. II VI The correlation functions used to characterize 
the dynamics are introduced in Sec.lVland finally the re- 
sults of the numerical calculations are discussed in Sec. 

ED 

II. HYPERFINE INTERACTION IN QUANTUM 
DOTS AND ITS SEMI-CLASSICAL 
REPRESENTATION 

The Hamiltonian describing the hyperfine coupling be- 
tween conduction band electrons and the lattice nuclei in 
GaAs is of the well known form of the contact potential 

H BF =AS-J2hS(r-R k ), (1) 

k 

where A is the hyperfine constant, S (Ik) and r (Rk) 
are respectively the spin and position of the electron (A:th 
nuclei). In the GaAs conduction band (which is mainly 
composed of s-orbitals), the dipole-dipole part of the hy- 
perfine interaction vanishes^. In this paper the focus 
will be on electrons localized in a quantum dot, and the 
hyperfine interaction in such systems. The quantum dots 
considered here are quite general, but we introduce some 
restrictions to simplify the model. First of all it is as- 
sumed that the number of electrons is fixed, preferably 
to one. From the experimental point of view this as- 
sumptions is quite reasonable since having only a sin- 
gle electron in the dot has already been demonstrated 
experimentally 2 ?'-. The second assumption is that the 
orbital level splitting is much larger than the hyperfine 
energy. In this case the hyperfine Hamiltonian can be 
projected to the lowest orbital level since contributions 



from higher orbitals are strongly suppressed due to the 
large orbital energy separation. If the ground state or- 
bital ip(r) of the quantum dot is known, then an effective 
spin Hamiltonian can be written as 

H s = gf x B B ■ S + 7 GaAs J2 B ■ tk + J2 A H(Rk)\ 2 S ■ I k 

k k 

(2) 

where g is the <?-factor, /is is the Bohr magneton, 7c a As 
is the gyromagnetic ratio of the effective nuclear species 
and B is the external applied field. 

In a typical quantum dot a single electron spin may be 
coupled to 10 4 — 10 6 nuclear spins. When the electrons 
interacts with so many nuclei it is possible to interpret 
the combined effect of the nuclei as an effective mag- 
netic field. Before proceeding further it is convenient to 
introduce a different way of writing the hyperfine inter- 
action in the last term in Eq. (J5J). The wave function of 
the ground state orbital has some characteristic spatial 
extent which is determined by the confining potential. 
Without loss of generality it may be assumed that the 
lateral (with respect to the underlying 2DEG) and trans- 
verse confining lengths are t and zq respectively. Defining 
the volume of the quantum dot as Vqd = ttzq£ 2 , a dimen- 
sionless function is introduced: 

f(R k ) = V QD \i;(Rk)\ 2 . (3) 

Furthermore, denoting the maximum value of / with 
/Max we introduce a dimensionless coupling constant 

g k =g(R k ) = l^le (0,1). (4) 

/Max 

The hyperfine coupling constant may be expressed in 
terms of the concentration, C„, of nuclei with spin / 
and a characteristic energy E n through the relation A — 
E n /C n I. The energy E n is the maximum Zccman split- 
ting possible due to a fully polarized nuclear system, its 
value being E n ps 0.135 meV in GaAs22i24i. The hyperfine 
interaction term can thus be written as 

H HF = S- K, (5) 

where we have introduced the operator for the effective 
nuclear magnetic field 

K = T^Sfclfc, (6) 

k 

and the characteristic hyperfine induced nuclear spin pre- 
cession frequency 



As was shown in Ref. it is possible to replace the op- 
erator in Eq. © with a effective nuclear magnetic field 
K . Its initial value is random since it is determined by 
unknown details of the nuclear system. Assuming that 
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the nuclear spin system temperature k B T 3> 7, all nu- 
clear states are equally likely to be occupied. Initially the 
nuclei are nominally decoupled from each other and the 
distribution of K is to a good approximation represented 
by a Gaussia n 12 ; 20 



P(K) = 



2it 1 2 NI 2 J 



.3/2 



exp 



3K 2 \ 



2<y 2 NI 2 J 



^ , (8) 



where I 2 = 1(1 + l)^ 1 £ fc g 2 . 

In the above discussion the dynamics of the nuclear 
spin system was disregarded. Before we include its dy- 
namics it is instructive to first derive exact operator equa- 
tions of motion, and then apply the semi classical approx- 
imation to those equations. The dynamics of the com- 
bined electron and nuclear spin systems are determined 
by the Heisenberg equation of motion 



d 

dt 
d 



S = K x S + gfi B B x S, 



(9) 



dt. 



Ik = IQkS X I k + 7GaAs-B X I k . (10) 



Multiplying Eq. (|10fl with gk and summing over k gives 
the equation of motion for K 



dt ' 



E 



9lh 



7GaAs# X K. 



(11) 



In contrast to the simple dynamics of Eq. [5J the equation 
of motion for K is quite complicated. The reason for the 
asymmetry is the position dependent coupling gk- The 
quantity in the brackets on the rhs of Eq. cannot 
be expressed in terms of K. Only in the simple case of 
constant gk it is possible to write a closed equation of 
motion for K and Without actually solving Eqs. 
@ and Ijlll) it is still possible to extract general fea- 
tures of the dynamics. For zero external magnetic field, 
the electron spin will precess with frequency cx E n N~ x l 2 
(which is the magnitude of the effective nuclear magnetic 
field) and the nuclear system precesses with frequency 
oc E n N~ x . Thus, for N 3> 1, the electron spin effectively 
sees a stationary nuclear system and in turn the nuclear 
system sees a time averaged electron spin. 

To incorporate (i) the separation of timescales and (ii) 
the inhomogeneous coupling we introduce a scheme that 
separates the nuclear system into subsystems, each 
being characterized by a fixed coupling g b . The effective 
nuclear magnetic field of a given subsystem is 



Kb = 19b ^2 



(12) 



keb 



where the notation k € b is shorthand for all nuclei whose 
coupling is g k € [gb - 5g/2,g b + Sg/2], with Sg = l/N b 
being the coupling constant increment. As long as N b <C 
N each subsystem can be replaced by a classical variable 
K b — > K b , which represents the effective nuclear field 
due to that particular nuclear spin subsystem. 



Using the same procedure as was used in deriving Eq. 
(|llf) we arrive at a equation of motion for K b and apply- 
ing the semi-classical approximation results in 



dK b 

dt 



jg b (S) x K b + 7caAs-B x K b , 



(13) 



where (S) is an appropriate time averaged electron spin. 
As will be discussed in the next section, this average elec- 
tron spin may be written as a function of the total effec- 
tive nuclear field 



K 



(14) 



The initial condition for each nuclear spin subsystem 
is randomly chosen from a Gaussian distribution whose 
variance is (see Appendix 



(Kl) 



7 2 iV/(7 + l)g b 8g. 



(15) 



The set of differential equations in Eq. (| 1 31) . with the 
associated random initial conditions constitute a set of 
autonomous differential equations. 

Separating the nuclear system into subsystem with a 
constant g b is an approximation to the continuous cou- 
pling gk- As the number of subsystems increases g b will 
more closely represent gk- However, for the semi-classical 
approximation to be valid each subsystem must contain 
many nuclear spins. Thus, increasing N b should better 
reproduce the actual system, as long as N b <C N. 



III. ADIABATIC APPROXIMATION FOR THE 
ELECTRON SPIN 

As we have shown in the previous section, the nuclear 
spin system may be treated as a slowly varying effective 
nuclear magnetic field acting on the electron spin. Let- 
ting H(t) represents any slowly varying magnetic field 
(fulfilling the usual adiabatic conditions) acting on a sin- 
gle electron spin, leads to the Hamiltonian 



H e (t) = 8 ■ H(t). 



(16) 



It is convenient to introduce the instantaneous eigenfunc- 
tions of the Hamiltonian, which are solutions of 



H e (t)\n(t)-±) =E±(t)\n(t);±) 



(17) 



The eigenstates are labeled by n(t) to indicate that these 
eigenstates are either pointing 'up' (+) or 'down' (-) 
along the total magnetic field, whose direction is deter- 
mined by the unit vector 



n(t) 



H(t) 
\H(t)\' 



(18) 



For a spin 1/2 in an external field, the eigenenergies are 
E±(t) — ±i|if(t)|, and the corresponding eigenstates 

are written in the basis of the S z eigenvectors 



!»(*);■ 
!»(*);■ 



1 



T> + o(*)li) 

^l + \a(t)\ 2 VU>- a *WIT> 



(19) 
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The time dependent mixing of spin components is given 
by 



o(t) = 



\H(t)\-H z (t) 
H x (t)-iH y (t)' 



(20) 



The wave function may be expanded in basis of instan- 
taneous eigenstates 



m)) = E c«(t)\n(t);a) (21) 

where the expansion coefficients are 

c±(t) =c±(io)e( l7±(t) -^ Jt o rfT£±(T) ). (22) 

The additional phase factor appearing in the previous 
equation is the usual adiabatic phased 

l±(t)=i [ dT{n{T):±\d/dT\n{T);±). (23) 



Integrating by parts the rhs of the last equation and using 
the orthogonality of the instantaneous eigenstates, it can 
be shown that the phase 7± is a real number. Also, for 
a doublet 7 = 7+ = — 7- • 

Using the wave function in Eq. I|21|l. the average elec- 
tron spin is 

(s(t)) = (m\s\m) 

= E \cAto)\ 2 (n(t);a\S\n(t);a) 

+ 25R| c ;(t )c_(t ) e ^ 2l7(t) " i ^ !iT|Jf(r)l )|24) 

The latter term in Eq. (|24H oscillates with frequency 
\H(t)\/h ^S> E n N^ 1 /h, so its average is zero on the 
timescales of the nuclear system. The average value of 
the electron spin entering Eq. (|13|) is 

(S) = ^2\ c M\ 2 (n(t)^\S\n(t);a) 

<T=± 

1 



(25) 



where 0q is the angle between the initial electron spin 
and n(0), see Fig. ^ The orientation between S and 
n(0) changes the precession of all K b 's by a multiplica- 
tive factor cos 9$. This overall factor has no effect on the 
dynamics and we subsequently put it to unity. Physi- 
cally, one might think of a dissipation mechanism that 
would initially align the electron spin and H to each 
other, although the mechanism itself is not critical for 
the following discussion. 

From the semi-classical version of Eq. (JHJ), the slowly 
varying magnetic field is H (t) = B + K(t) which results 
in an equation for (S) that depends only on K and B. 
Also, since we assume that the quantum dot is initially 
in the ground state orbital and that the orbital energy 
separation is much larger than the hyperfine energy, there 
are no 'Rabi oscillations' to higher orbitals. 



IV. SOME ANALYTIC PROPERTIES 

The rest of the paper will only consider the case of 
small magnetic field gfiBB <C K, which is the regime 
where the dynamics are most interesting. In the oppo- 
site situation the magnetic field strongly constraints all 
dynamics. The average electron spin, for B = is 



13)--- — 
{b >- 2\K\ 



(26) 



and the resulting equation of motions for the nuclear spin 
subsystems are 



— K b = jg b K x K b , 
at 



(27) 



where 7 = —j/2\K\, since \K\ is a constant of motion 



(28) 



and it equally changes the precession frequency of all the 
block K b . In addition to \K\, more integrals of motion 
can be constructed from Eq. I|27|l: 



= * w 

= d I= d( y K: 

dt dt \^ g b 




(29) 
(30) 

(31) 



The integral of motion in Eq. I|30|) . is actually the total 
spin I of the nuclear system- 2 ™. The integrals of motion 
are expected to affect the dynamics, i.e. the system will 
be 'constrained' by them. 

The solution for the electron dynamics is determined 
by the dynamics of the nuclear system, which is encapsu- 
lated in K[t). Although the dynamics are complicated 
there are some ways to characterize the motion of the 
nuclear system. For example, we can define the following 
quantity 



E 



1 



l+7bC 



K b (t) 



which satisfies the following equation of motion 
d „_ K(t;0)xK{t;Q 

dt [t ^'~ c 



(32) 



(33) 



It is possible to construct other similar equations, but in 
general no simple solution for them exist. 



V. CORRELATION FUNCTIONS 

A wide class of classical systems exhibits decaying cor- 
relation functions. This occurs in classically chaotic sys- 
tems in which the motion is such that the memory about 
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initial conditions is lost at some typical timescale^. Most 
of sufficiently complicated classical systems are eventu- 
ally chaotic. One might expect that the set of equations 
Eq. (|f 3[) should describe chaotic dynamics and decay of 
correlation functions. We will see later on that this is not 
the case. 

A useful way to characterize the electron spin dynam- 
ics it is to introduce certain correlation functions. For 
an isolated quantum system these correlation function 
are expected to oscillate periodically, without any decay. 
Incorporating environmental effects usually shows up in 
modified behavior of the correlation functions. The ex- 
pected behavior is that they should decay as a function 
of time. To investigate how the nuclear spin system acts 
as a spin bath (environment), we introduce the following 
correlation function 



G(t) = (t \S(t) ■ S\ t), 



(34) 



where the time evolution of the operators is in the usual 
Heisenberg picture. 

Since we are focusing on the slow dynamics it is useful 
to write these correlation functions for long timescales. 
In the adiabatic approximation the correlation functions 
may be written as 



G(t) = 



l (l-Kt)| 2 )(lH fl (0)| 2 )+4a(tK(0) 

4 (i + Ht)p)(i + | a (o)p) 



,(35) 



where the a's are defined in Eq. 1|2(J[1 . The most interest- 
ing regime corresponds to weak external magnetic fields. 
In this case there is no preferred direction and the dy- 
namics show the richest behavior. In that limit the cor- 
relation functions take the simplified form 



G(t) 



K(t) ■ K(0) 
4K(t)K(0) 



(36) 



From these equations it is evident that the electron spin 
correlation function is determined by the nuclear system 
variables for times t h\/N /E n . 

When dealing with many identical systems in which 
the electron can be prepared in a given initial state but 
the effective nuclear magnetic fields differ in the initial 
values, ensemble averaged correlation function must be 
considered. No information is available about the state 
of the effective nuclear magnetic field, except that their 
initial vales are Gaussian distributed. In this case the 
correlation function in Eq. I|35|) should be averaged over 
the appropriate distributions 



(G(t)) 



f l[dK bt0 P({K b>0 })G(t;{K bt0 }), 



(37) 



where K^.o = K b (0) and P is the Gaussian distribution 
of the initial values. Note that the correlation functions 
appearing in Eq. (|36|) are also functions of the set of 
initial conditions {K b (0)}. 



VI. RESULTS 

The correlation functions in Eqs. I|36[l an d (|37|l are 

in general not exactly solvable so numerical simulations 
have to be used. In order to calculate them time series 
for K b (t) need to be calculated. These are obtained by 
numerically integrating the differential equations in Eq. 
(I27[) . We focus on the case of no external magnetic field. 
In GaAs 7 « 10~ 7 meV~ 10 s Hz for quantum dots con- 
taining N w I0 6 nuclei. The differential equations are 
solved by integrating numerically Eq. (|27|) using the 4th 
order Runge Kutta method. The K b (t)'s are then used 
to calculate a(t) that enter Eq. (|35|) . 

For the ensemble averaged correlation functions many 
sets of time series {K b (t)} are calculated, each corre- 
sponding to different random initial condition chosen 
from a Gaussian distribution. The results for the en- 
semble averaged correlation function in Eq. I|37|l are pre- 
sented in Fig. |21 Each curve is the result of calculations 
for different number of subsystems N b = 4, 8, ... , 512. 
As is to be expected, the correlation functions decay in 
time but a saturation value is reached for sufficiently long 
times, which is determined by t sat oc 7 _1 N b . This satura- 
tion is an artifact of the discretization, i.e. it introduces 
a time above which the calculated correlation function 
no longer represents the true correlation function. For 
the calculated correlation function to have a meaningful 
limit, the saturation value should approach zero as N b in- 
creases and the saturation time should go to infinity. The 
inset in Fig. [21plots the saturation values of (G(t)). It is 
evident that the saturation values tend to zero for larger 
N b . The decay fits an inverse logarithm a/ \n(f3N b ) quite 
well. This indicates that there is a well defined N b — > 00 
correlation function which has an inverse logarithmic de- 
cay oc 1/lni as i — > 00. The origin of the logarithmic 
decay is not fully understood (and thus the theoretical 
values of a and (3) but it might be related to the total 
spin \I\ 2 oc lni sat (see Appendix 1X1) . The correlation 
function is normalized and its value might be dominated 
by 1/|/| 2 fort>i sat . 

The correlation functions for a single system, i.e. with- 
out taking the ensemble average, yield quite different re- 
sults. The calculations are performed for N b = 8,32, 128 
and 256. For each number of subsystems N b the calcula- 
tions were repeated for various random initial conditions, 
but no averaging is performed. The results of the calcu- 
lations for G(t) are presented in Figs. 13161 (note different 
range on the horizontal, or t-axis). The common feature 
of all the curves, for all values of Nb, is that they do not 
decay with time. This behavior persists to even longer 
times, not shown here. Even though more complicated 
behavior is observed for large N b , the characteristic fre- 
quencies of the correlation function oscillations do not 
show any obvious dependence on the number of subsys- 
tems. 

It is instructive to look at the power spectrum of G(t), 
where the oscillating behavior of the correlation func- 
tions becomes more apparent. The power spectrum (or 
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the squared Fast Fourier Transform) of G(t) is shown for 
the 1st and 6th curves (counted from the bottom one) 
for N b = 32 (see Fig. EJ), in Figs. and El respectively. 
The sharp, isolated peaks in the spectra illustrate well 
the multi-periodic oscillations observed in G(t). This be- 
havior is still present in the power spectra for N b = 256. 
Figs. l9llTl corresponding to respectively the 2nd, 4th and 
6th curves in Fig. [5] show that even for such a compli- 
cated system (256 coupled, non-linear differential equa- 
tions) the correlation functions still show sharp, isolated 
peaks corresponding to well defined oscillation periods 
and additional many smaller, closely spaced peaks. 

The power spectrum is the Fourier spectrum of the 
time averaged correlation functions, which are completely 
different from the power spectra expected for the ensem- 
ble averaged correlation functions shown in Fig. J2J). This 
implies that the time average and the ensemble average 
are not equivalent, i.e. the system in question is not er- 
godic. The simplest way to think about this is to con- 
sider the integrals of motion for the system. In the case 
of the time-averaging the motion of the system is at all 
times 'constrained' by the integrals of motion, resulting 
in multi-periodic correlation functions that show no de- 
cay in time^. For the ensemble case, the averaged cor- 
relation function get contributions from many 'systems' 
which have different values of the integrals of motion that 
results in an effective cancellation of periodic oscillations, 
leading to a decay of the correlation functions. 

The behavior of the correlation functions for the single 
system may be broadly explained in this way: The oscil- 
lations of the correlation function reflect that the system 
is in some sense close to being exactly solvable. These 
features will probably vanish if further terms are included 
into the Hamiltonian in Eq. (J2J). The most natural term 
would be the dipole-dipole interaction between the nu- 
clei which would kill most of the integrals of motion. It 
is important to recognize that the timescale related to 
the dipole-dipole interaction is very long, of the order 
l(U 3 s, and the resulting decay time would reflect that. 

In connection to coherently controlling the spin, the 
motion of the electron spin in the effective nuclear mag- 
netic field will cause 'errors'. Even though in this model 
nuclear spins do not decohere the electron spin (in the 
sense of not causing decay of correlation functions as a 
function of time), they leads to a complicated, and un- 
predictable, evolution. Consequently, an electron spin 
initially in | f ) can be found in the opposite spin state 
on a timescale oc Thus, even though the hyperfme 

coupling to the nuclear system does not lead to decoher- 
ence in this model it can strongly affect the dynamics of 
the electron spin. 

The authors acknowledge financial support from FOM 
and SIE would like to thank Oleg Jouravlev, Dmitri 
Bagrets and Lieven Vandersypen for fruitful discussions. 



APPENDIX A: THE VARIANCE IN TERMS OF g b 

The initial condition for each nuclear spin subsystem 
is chosen from a Gaussian distribution whose variance is 
determined by 



(K 2 ) = 7 2 I(I+l)g 2 C n V b . 



(Al) 



where V b is the volume of subsystem &, and C n V b is the 
associated number of nuclear spins in that volume. For 
the calculations it is convenient to express the variance 
for a given subsystem only in terms of the coupling g b by 
expressing the subsystem volume V b as a function of g b . 
The volume of the subsystem V b is related to g b via 



dV 
dg 



89, 



(A2) 



where V is the volume of the region where g > g b . 
The functional form of g(r) is determined by the density 
|^(r)| 2 . In the numerical calculations a lateral parabolic 
confinement of an underlying 2DEG is used. Assuming 
a constant electron density in the z direction (growth 
direction) the electron density is 



g(r) = 9(z /2-\z\)eM-(x 2 + y 2 )/t 2 ), 



(A3) 



which gives the simple relation x 2 + y 2 = £ 2 \n(l/g), 
within the 2DEG. Using the relation for the volume 
V(r) = zqtt(x 2 + y 2 ) the subsystem volume is 



Vb = U QD ^, 
9b 



(A4) 



which gives the variance of the distribution of the effec- 
tive nuclear magnetic field for a given subsystem in terms 
of the coupling 



(Kl) 



j 2 NI{I + l)g b 5g. 



(A5) 



These results can be used to calculate the variance of 
the total nuclear spin defined in Eq. I|3U|I 



(Kb ■ K w 

9b9v 



E 

b.V 

(K) ^ Sg 
9 2 ^9b 



(A6) 
(A7) 



« ^ln(2JV 6 ), 7V b »l. (A8) 
9 

In the last step it is assumed that g b = 1 — (b — 1/2) /N b , 
resulting in the logarithm. 



* Present address: Department of Physics and Astron- Switzerland 
omy, University of Basel, Klingelbergstrasse 82, CH-4056, 
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FIG. 1: The time dependent electron spin (S(t)) precesses 
rapidly around the total effective magnetic field, resulting in 
a slowly varying average spin {S) that the nuclei see. The 
angle between the instantaneous electron spin and n is do 
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FIG. 2: The ensemble averaged correlation function as a 
function of time for N b = 4,8,16,32,64,128,256 and 512. 
The inset show the asymptotic values of (G(t)) and a fit to 
a/\n(f3N b ). 
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FIG. 3: Numerical calculations of the correlation function 
G(t) for Nb — 8 and various randomly chose initial conditions. 
The curves are offset for clarity and the vertical range is the 
same for all curves, i.e. —0.25 to 0.25 
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FIG. 4: Numerical calculations of the correlation function 
G(t) for Nb = 32 and various randomly chose initial condi- 
tions. The curves are offset for clarity and the vertical range 
is the same for all curves, i.e. —0.25 to 0.25 
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FIG. 5: Numerical calculations of the correlation function 
G(t) for Nb — 128 and various randomly chose initial condi- 
tions. The curves are offset for clarity and the vertical range 
is the same for all curves, i.e. —0.25 to 0.25 
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FIG. 6: Numerical calculations of the correlation function 
G(t) for Nt — 256 and various randomly chose initial condi- 
tions. The curves are offset for clarity and the vertical range 
is the same for all curves, i.e. —0.25 to 0.25 
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FIG. 7: The power spectrum of the 2nd curve (counted from 
the bottom one) in Fig. [I] corresponding to Nb = 32. 
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FIG. 8: The power spectrum of the 6th curve (counted from 
the bottom one) in Fig. [I] corresponding to Nb = 32. 
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FIG. 9: The power spectrum of the 2nd curve (counted from 
the bottom one) in Fig. corresponding to Nb = 256. 
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FIG. 10: The power spectrum of the fourth lowest curve in 
Fig. corresponding to Nb = 256 
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FIG. 11: The power spectrum of the top curve in Fi 
corresponding to Nt = 256 



